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Evaporative cooling of helium nanodroplets is studied with a statistical rate model that includes, 
for the first time, angular momentum conservation as a constraint on the accessible droplet states. 
It is found that while the final temperature of the droplets is almost identical to that previously 
predicted and later observed, the distribution of total droplet energy and angular momentum states 
is vastly more excited than a canonical distribution at the same temperature. It is found that the 
final angular momentum of the droplets is highly correlated with the initial direction, and that a 
significant fraction of the alignment of the total angular momentum should be transferred to the 
rotational angular momentum of an embedded molecule. 



The study of helium nanodroplets has received consid- 
erable attention in the past decadeiS^ Such droplets 
rapidly cool by helium atom evaporation while they 
travel inside a vacuum chamber. Brink and Stringar: 4 
used a statistical evaporation model and predicted a ter- 
minal temperature of He nanodroplets of 0.4 K, in ex- 
cellent agreement with the value of 0.38 K later deduced 
from the rotational structure in vibrational spectra of 
SF 6 and other molecules^ Despite the obvious success 
of this theoretical work, the model used is clearly incom- 
plete in that the constraint of angular momentum conser- 
vation was not imposed. The need for a more complete 
evaporative cooling study was made evident by the recent 
observation of a polarization anisotropy in the absorption 
spectrum of pentacene in helium droplets^ The authors 
of this work suggested that the total angular momen- 
tum deposited in the droplets by the pickup of a heavy 
molecule is aligned perpendicular to the droplet velocity, 
and that this droplet alignment survives the evaporative 
cooling and is transferred to the embedded molecule. The 
present study was undertaken to test the reasonableness 
of this conjecture. 

We model the evaporative cooling with a statistical 
rate approach, analogous to phase space theory in uni- 
molecular dissociation, which explicitly includes the con- 
straints of angular momentum conservation^ We use 
Monte Carlo sampling to follow cooling 'trajectories' as 
the droplets lose much of their initial energy and angu- 
lar momentum by helium atom evaporation. It is found 
that the droplets cool to final temperatures close to those 
predicted without angular momentum constraints. 4 How- 
ever, the distribution of terminal droplet states (where 
the evaporative cooling lifetime becomes longer than typ- 
ical flight times in experiments) cover a vastly broader 
range of energy (E) and total angular momentum (J) 
than was previously expected. Further, it is found that 
the final angular momentum vector of the droplet is 
highly correlated with the initial value, such that much 
of the alignment remains, and that a sizable fraction of 
this alignment is transfered to an embedded rotor. 

Evaporative Cooling Model — Consider a helium nan- 
odroplet, D, with initial values of the conserved quanti- 
ties n (number of helium atoms), E' (total internal en- 
ergy in units of Kelvin), and J' (the total angular mo- 



mentum in units of h). If E' is sufficiently great, the 
droplet will cool by helium atom evaporation by the re- 
action: D(n, E', J') -> D(n - 1, E", J") + He (^trans, L) 
where -Etrans is the center of mass translational energy in 
the dissociation and L is the orbital angular momentum 
quantum number of the fragments. Conservation of total 
angular momentum requires that J', J" , L obey the tri- 
angle rule. Based upon the bulk density of helium^ the 
droplet has radius R(n) — r^n 1 ^ with ro = 2.22 A. Con- 
servation of energy, including the change in the surface 
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where E^ = 7.2 K is the binding energy of helium to 
the bulkiS Because of the centrifugal barrier, classically 
£tra„s > f? L{L + 1)/ (2m He R(n) 2 ) , leading to a E^ ax 
value for each L. 

Using statistical reaction rate theory, the rate of he- 
lium atom evaporation is given by& 
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p(n, E', J') is the density of states of the droplet at fixed 
angular momentum, and N a (n, E' , J') equals the number 
of states of the product droplet and departing He atom 
consistent with total E' and J' and with the departing He 
atom having kinetic energy above the centrifugal barrier. 
If we denote the total number of n atom droplet states 
with energy < E and total angular momentum quantum 
number J by N(n, E, J), we can use the triangle rule to 
write N {n,E',J') as a sum of N(n - 1, E^ ax (L), J") 
over all allowed values of L,J". 

For droplets of the size and energy range of interest 
to experiments, the density of states is dominated by 
quantized capillary waves on the surface of the droplets 
( ripplon modes). 4 Simple, but highly accurate, analyt- 
ical approximations for N(n, E' , J') and p(n,E',J') = 
(dN/dE)j, have recently been published^ thus allow- 
ing calculation of k(n, E, J). Using an ensemble that 
conserves n, E, J, the droplet temperature can be cal- 
culated using (ksT(n, E, J)) -1 = (<9lnp(n, E, J)/dE)j. 
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For droplets with only ripplon excitations, the densities 
of states depends upon n only through the size of a re- 
duced energy scale which is equal to 3.77n _1 / 2 K. Start- 
ing with a given initial size, energy, and total angular 
momentum for the droplet at time zero, we calculate 
the evaporation rate and advance time by the inverse 
of this number. We then randomly pick a single open 
decay channel (E",J",L) with probability proportional 
to N(n — 1, E", J"). We can treat the product droplet as 
a new initial condition, and calculate another evapora- 
tion event. As the droplet cools, the rate of evaporation 
falls exponentially. When the cumulative time for the 
next evaporation event is greater than the assumed free 
flight time in an experiment, we terminate the evapora- 
tion process. 

We also did evaporative cooling calculations for 
droplets with an embedded linear rotor with effective 
rotational constant B, which will typically be several 
times smaller than the corresponding value in the gas 
phase due to helium motion induced by the rotation^ 
For these calculations, the integrated density of states, 
N c (n 7 E,J), for the combined droplet + rotor system 
is calculated by convolution of the integrated droplet 
density of states with the spectrum of the rigid ro- 
tor: N c (n, E, J) = £ . £jh N(n, E - Bj(j + 1), J') 
These were calculated for various values of B and the 
resulting densities fit to the same functional form pre- 
viously used for N(n,E,L), but with the fitting con- 
stants now expressed as polynomial function of B and 
its inverse. Differentiation gives the combined density of 
states, p c {n, E, J). 

Alignment is defined as (P2(cos0)), where 9 is the an- 
gle of the angular momentum vector with the quantiza- 
tion axis (chosen to be the velocity vector of the helium 
droplet beam) . In order to track the changes in alignment 
during cooling, we must propagate the angular momen- 
tum projection quantum numbers of droplet (M) and 
embedded rotor (to). The final alignments are found to 
be proportional to the initial alignment of the total an- 
gular momentum created by the pickup process. 

If we have an initial distribution P(n, M') over 
angular momentum projection quantum numbers M' , 
we can use the Wigner-Eckart theorem to cal- 
culate the product distribution of the projection 
quantum number populations^ P(n — 1,M") — 

J2m> ( J "> M "> L > M> ~ M"\J', M') 2 P(n, M'). The as- 
sumption of equal probability of all states consistent with 
total E and J leads to a similar expression for the prob- 
ability of populating a rotor quantum state j, m. This 
allows the calculation of the alignment of the rotor an- 
gular momentum, which is an experimental observable. 

Results — We did evaporative cooling calculations with 
initial values of n = 10 4 , E = 1700. K, and J = 
10, 1000,2000, ...5000, computing 2500 Monte Carlo 
'trajectories' in each case. These values where selected as 
likely initial conditions following the pickup of a tetracene 
moleculei 1 ^ These conditions correspond to an initial 
droplet temperature of 2.61 K independent of J. We also 



assumed that the initial condition for the total angular 
momentum is P(M) — Sq,m (initial alignment = —1/2). 
A cooling time of 100 (is was used. 

We can first make some qualitative predictions. We 
expect that on average the evaporated helium atoms 
will reduce droplet energy by E^ + (3/2)k&T where T 
is the average temperature during evaporation, which 
would imply an energy loss of « 10 K per He evapo- 
rated. Thus, we would expect ~ 170 helium atoms to 
be evaporated. A helium atom evaporated at T will 
carry root mean squared (rms) angular momentum of 
y/2mH e kBTR ~ 23h. If the departing atom's angular 
momentum is parallel to the initial angular momentum 
(as for water droplets leaving a wet, spinning tennis ball), 
we can expect the droplet to lose about 4000?i while 
cooling. Alternatively, if the orbital angular momenta 
of the evaporated atoms are random, the droplet angu- 
lar momentum would undergo a 3D random walk, and 
the mean squared angular momentum of the droplets 
would be expected to increase by the number of evap- 
orated atoms times the square of the angular momentum 
per evaporation event. The form for p(n, E, J) for fixed 
n, E is that of a thermal distribution of a spherical rotor 
with a mean value of J( J + 1) that grows as E 8 /"' iii For 
the initial values of n, E, the rms of the distribution of 
p(n, E, J) occurs for J = 345, an order of magnitude be- 
low the mean value of the angular momentum deposited 
in pickup of a large molecule, such as pentacene. It can 
be expected that during evaporation, the density of states 
will bias evaporation such that the J moves towards the 
rms value, which will result in the average angular mo- 
mentum being reduced during evaporation. For a canon- 
ical distribution at 0.38 K, the mean terminal value of 
(E) = 18.5 K and y/(J(J+ 1)) = 39 are predicted for 
droplets of n = 10 4 iii 

Figure ^ shows a plot of droplet energy, angular mo- 
mentum quantum number, temperature, and number of 
evaporated helium atoms as a function of time for five 
representative trajectories with initial condition Ei nt = 
1700. K and J int = 3000. It is found that most of 
the evaporation takes place in the first ~ 1 ns, but 
that the temperature slowly drops even for long times. 
The final temperatures of each trajectory is similar to 
that observed experimentally for molecules inside he- 
lium droplets*^ and also to that predicted by evapo- 
rative cooling calculations without angular momentum 
constraints^ The remarkable difference from those pre- 
vious calculations is that the final energies and angular 
momenta are vastly higher than that found in the previ- 
ous work. 

Table I gives a summary of the numerical results for a 
range of initial angular momentum values. For calcula- 
tions with a shorter cooling time (10 fis), the results were 
almost identical, except that the final temperatures were 
systematically w 50 mK higher. It is seen that droplet 
cooling gives a mean residual energy and angular mo- 
mentum vastly higher than that of a canonical ensemble 
at the same final temperature, and this energy and angu- 



3 



E 



3500 
3000 
2500 
2000 
1500 
1000 
500 

V 

















5 < ^ 






































D D 








#838 8 






































k ,A . 














1 1 1 


□nib mi 



10 11 Iff 5 KT 7 TO 5 

Evaporation Time / sec 



175 




150 










£ 


125 


o 




ns 




<D 


100 


I 




"D 




.2 


75 


2 




o 




Q. 


50 


CO 

> 

LU 




* 


25 










FIG. 1: Five evaporative cooling trajectories starting with 
initial conditions J = 3000 ft, E = 1700. K, n = 10 4 , which 
predicts an initial droplet temperature of ~ 2.5 K. Plotted 
are the droplet angular momentum, temperature, energy, and 
number of helium atoms evaporated as a function of time since 
the start of the simulation 



lar momentum rises quickly with initial angular momen- 
tum value. Notice that the 'trapping' of angular momen- 
tum is found even for initial values considerably smaller 
than the 'maximum' value predicted above that can be 
shed by evaporation. It is also evident from the table 
that most of the initial alignment of the total angular 
momentum is maintained which indicates that the final 
angular momentum is nearly parallel to the initial angu- 
lar momentum. While for each initial condition, a broad 
distribution of final E and L values are found, these are 
distributed in a narrow band of energy width w8K, fol- 
lowing a line corresponding to constant temperature. For 
fixed J, the final E values are ss 45% larger than the low- 
est possible droplet energy for that J, which corresponds 
to J/2 quanta in the lowest, L = 2 ripplon mode. 

While a systematic study of the final distributions 
upon the full range of likely droplet initial conditions 
is beyond the scope of the present report, we would like 
to indicate the trends. If the initial energy and total an- 
gular momentum quantum number are kept constant at 
1700. K and 1000 respectively, and the initial size of the 
droplet is decreased from n — 10 4 to 3 x 10 3 , the initial 
temperature rises from 2.61 to 3.69 K. As expected, the 
higher initial temperature rises the average kinetic energy 
carried away by the He atoms (3.53 — > 4.88 K), resulting 
in a modest decrease (164 — > 149) in the average number 
of evaporated He atoms. Because of decreased droplet 
size, the rms orbital angular momentum carried away by 
the He atoms decreases (25.7 — > 20.4). Despite this, the 
final average angular momentum decreases (401 — > 355), 
due to an increased correlation in the emission direction 
(average cosine of angle between L and J increases from 
0.156 to 0.23). The final average energies and tempera- 
tures are almost unchanged. If we keep the initial angular 
momentum at J = 1000 and the size at n — 10 4 but de- 
crease the initial energy from E — 1700. to 1000. K, we 
find that the initial temperature (2.61 — > 2.08 K), and 



number of evaporated helium atoms (164 — > 101) both 
decrease as expected. The He atoms mean kinetic en- 
ergy (3.53 — > 2.87 K) and rms orbital angular momentum 
(25.7 — > 23.0) both decrease as expected for lower initial 
temperature. Unexpectedly, the average final droplet en- 
ergy (44.2 — ► 49.8) and angular momentum (401 — > 475) 
increase when the initial energy is decreased. 

We performed evaporative cooling calculations for B = 
0.6 and 1.2 GHz. These values were selected as they rep- 
resent typical values for "heavy rotors" that have been 
studied in helium nanodroplets with rotational resolu- 
tion^ As expected, the presence of the rotor has lit- 
tle effect upon the final distribution of E and J of the 
droplets. Figure |2 shows the calculated ratio of align- 
ment of the rotor angular momentum to that of the total 
angular momentum for droplets with the average final 
values of E and J for initial states corresponding to the 
J initial = 10, 1000, 2000 . . .5000. It can be seen that the 
degree of rotor alignment increases with both j and the 
initial angular momentum, though the alignment ratio 
appears to saturate at higher Jinitial values. The level 
of alignment found in this work is certainly within the 
range that should be detectable by experiments of the 
type reported for pentacene (« 10 — 20%). 7 
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FIG. 2: Plot of the ratio of the alignment of an embedded 
rotor to the alignment of the total angular momentum of the 
droplet, as a function of the rotational quantum number of 
the rotor. Curves are for final states reached from initial 
conditions E = 1700. K and J = 1000, . . . 5000 starting from 
the lowest curve 

Experimental Consequences- This work has demon- 
strated that the distribution of the internal excitations 
of helium nanodroplets should be vastly more excited 
than had been previously predicted, based solely upon 
the low temperature of the droplets. The present re- 
sults support the interpretation of Portner et at, that a 
remnant of the initial pick up angular momentum sur- 
vives evaporative cooling and provides a bias that par- 
tially aligns embedded molecules i The present work pre- 
dicts that this should be a common phenomenon. It also 
provides a rationalization for the failure of a previous at- 
tempt to predict the rovibrational lineshapes of embed- 
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TABLE I: Evaporative Cooling of Helium Nanodroplets initially with 10000 helium atoms and initial internal energy equal to 
1700. K. "±" after each entry gives the standard deviation of the ensemble of 2500 cooling trajectories. 



Jinitial 


10 


1000 


2000 


3000 


4000 


5000 



(An) 166.05 ± 3. 15 164.39 ±3. 09 160. 75 ± 2. 86 157.24 ± 2. 83 153.63 ± 2.68 150.05 ± 2.46 

(Esnex) 25.88 ±5.66 44.45 ± 8.56 81.06 ±9.96 125.0 ± 10.1 174.98 ± 10.30 229.2 ± 9.9 K 

(J fina i) 170. ± 64. 404 ± 107. 897. ± 135. 1528. ± 145. 2273. ± 152. 3108. ± 151. 

(T fina i) 354.7 ±21.9 347.1 ± 17.2 343.4 ± 11.4 342.8 ± 8.2 344.4 ±6.2 346.1 ±5.0 mK 

(J Sttags) 3.53 ± 2.79 3.52 ±2.75 3.53 ± 2.71 3.47 ± 2.65 3.39 ± 2.56 3.26 ± 2.44 K 

y/{L(L + l)) 25.6 ± 26.2 25.7 ± 26.0 26.1 ± 26.1 26.5 ± 26.1 26.6 ± 25.9 26.6 ± 25.6 

(cos(6»j/ jL )) 0.03 ±0.58 0.156 ± 0.572 0.267 ± 0.554 0.357 ± 0.530 0.426 ± 0.504 0.480 ±0.480 

{ /(final) alignment) q.002 ± 0.020 0.829 ± 0.053 0.958 ± 0.007 0.983 ± 0.002 0.992 ± 0.001 0.995 ± 0.0004 

J (initial) alignment 



ded molecules based upon the inhomogeneous distribu- 
tion of molecule "particle in a box" statesJ^ That work 
assumed that these translational states followed a Boltz- 
mann distribution, while the present work suggests that 
high angular momentum translational states of the em- 
bedded molecules will be substantially more populated 
in the droplets. More significantly, the dramatic effect of 
angular momentum constraints calls into question the use 
of Path Integral Monte Carlo simulations for nanodroplet 
experiments since such calculations assume a canonical 
ensemble. The present work demonstrates droplets pre- 
pared in existing experiments are poorly described as a 
canonical ensemble. 

Many of the predictions of the present calculations 
could be directly experimentally tested. The trapping 
of increasing energy in the droplets with increasing ini- 
tial angular momentum will reduce the expected amount 
of helium atom evaporation as the collision energy is 
changed. However, it is not clear if this could be dis- 
entangled from other effects, including a collision and 
impact parameter dependence of the pickup probabil- 
ity and the possibility of coherent ejection of helium 
atoms during the pickup process, i.e. the early evapo- 
ration events could be strongly nonstatistical and thus 
poorly predicted by the present model. Measurement 
of the spatial distribution of helium atoms evaporated 



from a droplet would be very revealing, as the present 
calculations predict that the atomic velocity distribution 
is highly anisotropic in the center of mass frame. Since 
the orbital angular momentum of the evaporated atoms 
is on average parallel to the initial angular momentum, 
which itself is largely perpendicular to the droplet ve- 
locity, one expects to find the evaporated helium atoms 
peaked ahead and behind the droplets in the center of 
mass frame. Measurement of the radial distribution of 
atoms or molecules embedded in helium nanodroplets 
would also reveal the expected increased translational 
angular momentum of the dopant. It is possible that 
high energy electron or X-Ray scattering experiments 
could reveal the radial distribution function of atoms or 
molecules solvated by helium nanodroplets. The most 
direct test of the present model would be to measure the 
trapped internal energy and/or angular momentum of 
the droplets. One way that could be done in principle is 
to measure the Stokes and antistokes Raman spectrum 
from the lowest, L — 2 ripplon mode. It should be noted 
that excitation of the L — 2 mode should dominate in 
the trapped states, which have nearly the lowest possible 
internal energy for a given total angular momentum. 
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